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Abstract: The effect of a temperature dependent bulk viscosity to entropy density ratio (C/ s ) along with a constant 

shear viscosity to entropy density ratio {rj/s) on the space time evolution of the fluid produced in high 
energy heavy ion collisions have been studied in a relativistic viscous hydrodynamics model. The boost 
invariant Israel-Stewart theory of causal relativistic viscous hydrodynamics is used to simulate the evolution 
of the fluid in 2 spatial and 1 temporal dimension. The dissipative correction to the freezeout distribution 
for bulk viscosity is calculated using Grad's fourteen moment method. From our simulation we show that 
the method is applicable only for (/s < 0.004. 

PACS (2008): 12.38.Mh ,47.75. +f, 25.75.Ld 

Keywords: Bulk Viscosity • Relativistic hydrodynamics • Grad's moment 
© Versita Warsaw and Springer- Verlag Berlin Heidelberg. 



1. Introduction 

Recent experiments in high energy nuclear collisions at relativistic heavy ion collider (RHIC) confirms the existence 
of a new state of matter known as Quark Gluon Plasma (QGP) [1]. The production of QGP in heavy ion collision 
and its subsequent collective evolution provide us the unique opportunity to study the transport properties of this 
most fundamental form of matter. Relativistic viscous hydrodynamics simulations of observables like elliptic flow 
(^2) and transverse momentum (pr) spectra have been compared to experimental data to extract the QGP rj/s. 
Most studies show that the estimated value of rj/s lies between 1 — 4 x (1/4tt). However to correctly extract the 
rj/s of the QGP fluid, it is important to know the effect of finite bulk viscosity on the fluid evolution. Theoretical 
calculations based on pQCD [2] and lattice QCD [3] shows that the bulk viscosity is non-zero for the temperature 
range applicable in the heavy ion collision. In this work we use a temperature dependent form of C/s to study the 
effect of bulk viscosity in fluid evolution. The dissipative correction to the freezeout distribution function bulk 
viscosity has also been considered using Grad's 14 moment method. 
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Figure 1. (Color online) Q/s as a function of temperature. Red dashed line is the 77/s = 1/47T. 

2. Viscous hydro dynamic model 

The space time evolution of the fluid was simulated by simultaneously solving the energy momentum conservation 
equation d^T^ — 0, along with the relaxation equation for shear and bulk stress. Here T^ v — (e + _p + W)u^u v — 
(p + H)g^ u + tt^ u is the energy momentum tensor and e,p are energy density, pressure of the fluid, is the 
metric tensor; II and ty^ u are bulk and shear stress tensor respectively. According to the Israel- Stewart theory of 
causal viscous hydrodynamics [4], the shear and bulk viscosity obey the following relaxation equations Dty^ = 



^[2ryV <M "" > - ttH - (u^tt uX + u u ^ x ) Du X ] and DU = ~^[U + CV >^ + §C™9 M (^fr^)]. Here D is the 



convective derivative,^ and rn are the relaxation time for shear and bulk stresses respectively. We assume that 
the fluid achieve near local thermalization at proper time 0.6 fm. Initial transverse velocity (vt) is assumed to 
be zero. The initial energy density profile in transverse plane is calculated from a two component Glauber model 
with a central energy density eo = 30GeV/ fm 3 . Initial value of tt^ and n was set to their corresponding Navier- 
Stokes estimate. We assume the fluid freezes out when an element of it cools down below a constant temperature 
Tf = 130 MeV. The freezeout procedure was carried out by using Cooper-Frey algorithm. In the present study, 
we have used an equation of state (EoS) where the Wuppertal-Budapest lattice calculation [3] for the deconfined 
phase is smoothly joined at crossover temperature 174 MeV, with hadronic resonance gas EoS comprising all the 
resonances below mass ra re s=2.5 GeV. Q/s and 77/s are inputs to viscous hydrodynamics simulation. Figure 1 
shows the C/s(T), where (/s in the QGP phase is obtained by using pQCD formula (/s = 15 J (T) (1/3 - c 2 s (T)) 2 , 
the squared speed of sound c 2 s was calculated from lattice data [5]. In the hadronic phase Q/s is parametrized 
from [6]. The red dashed line in figure 1 is 77/s. 



We first discuss the change in pion pr spectra and V2 due to bulk and shear viscosity in the fluid evolution only. 



3. 



Results and discussion 



In the left panel of figure 2 temporal evolution of spatially averaged transverse velocity ((vt)) = 



is 



((7)) 
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Fi gure 2. (Color online) The left plot shows the temporal evolution of fluid transverse velocity for ideal (red), bulk (dashed dot) 
and shear viscous (dotted) evolution. The right plot is the 7r _ invariant yield as a function of pt for ideal (red) 
,bulk dashed dot and shear viscous (dotted) evolution. The inset figure shows the relative correction to invariant yield 
due to the bulk viscosity in comparison to ideal fluid. 
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Fi gure 3. (Color online) The left plot is the temporal evolution of momentum anisotropy of the fluid for ideal(red),bulk(dashed 
dot) and shear viscous(dotted) evolution. Right plot is the corresponding elliptic flow (i^) of tt~ . The inset plot shows 
the relative change in V2 due to bulk viscosity compared to ideal fluid. 



shown for ideal, shear and bulk viscous fluid. Here the angular bracket denotes space average and 7 = > 

\/l — v 2 

Because of the reduced pressure in bulk viscous evolution, {(vt)) is reduced in comparison to ideal fluid evolution. 
Whereas shear viscosity increase the pressure in the transverse direction, as a result the {(vt)) is larger compared 
to ideal fluid. The effect of the changed fluid velocity in viscous evolution is reflected in the slope of the pr 
spectra of tt~ shown in the right panel of figure 2. The relative change in the tt~ invariant yield (SN/Nideai, SN = 
Nbuik — N ideal) due to the bulk viscosity in comparison to ideal fluid is shown in the inset of right plot of figure 2. 
The relative correction is within ~ 10%. 

The temporal evolution of momentum space anisotropy e p = j a^%^ xx +ryy] * s shown in the left panel of figure 3. 
Viscosity tries to diminish any velocity gradient present in the fluid, as a result of that e p is smaller for both shear 
and bulk viscous evolution compared to ideal fluid. In a hydrodynamic model V2 is proportional to e p hence a 
reduction in e p will result in a reduced V2> Elliptic flow of tt~ as a function of pr is shown for ideal, shear and 
bulk viscous evolution in the right plot of figure 3. The inset shows the relative correction to V2 due to bulk 
viscosity. The relative correction to V2 is within ~ 3%. 
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Figure 4. (Color online) The invariant yield of tt~ as a function of pt for ideal and bulk viscous evolution with four different Q/s 
values are shown in the upper left plot. The lower left plot shows the relative correction to the tt~ yield as a function 
of pt for four different C/ s - See text for details. The right side plot is same as the left but for v<i. 

We have employed Grad's fourteen- moment method for calculating the dissipative correction to the freezeout 
distribution function as described in [7]. The details of the implementation of this method to our viscous 
code " ' AZHYDRO-KOLKATA" ' can be found in [8, 9]. The top left panel of figure 4 shows the p T spectra 
of pions for ideal(black solid line) and for four different values of £/s. The corresponding relative correction to 
the pt spectra is shown in the bottom left panel. The V2 of pion and the relative correction is shown in the right 
panel of figure 4. Freeze-out correction in Grad's moment method is obtained under the assumption that the 
non-equilibrium correction to the distribution function is small than the equilibrium distribution function. It is 
then implied that the relative correction SN/N eq is small for Grad's method to be applicable. The shaded band in 
the bottom left panel corresponds to the relative correction of 50%. If we consider here a correction of magnitude 
greater than 50% indicates the breakdown of the the freezeout correction procedure then our study shows that 
the Grad's method will be applicable if the (,/s has value less than 0.01 times the present form considered here. 
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